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Abstract 

The beam-beam effects of the PEP-II as an asymmetric 
collider are studied with strong-strong simulations using a 
newly developed particle-in-cell (PIC) code[[l|]. The simu- 
lated luminosity agrees with the measured one within 10% 
in a large range of the beam currents. The spectra of co- 
herent dipole oscillation are simulated with and without 
the transparency symmetry. The simulated tune shift of 
the coherent tt mode agrees with the linearized Vlasov the- 
ory even at large beam-beam parameters. The Poicare map 
of coherent dipole is used to identify the beam-beam reso- 
nances. 

1 INTRODUCTION 

The PEP-II |j| and KEKB [§] as asymmetric collider, 
which consists of two different rings at different energy, 
have been successfully constructed and fully operational. 
The beam-beam effects in this new type of e + e~ collider is 
one of the important physical phenomena to be studied be- 
cause, with twice of more parameters, there are much more 
choices of operating parameters to gain a higher luminos- 
ity. Basically, there are two major choices for the operat- 
ing parameters. One choice is the symmetric parameters 
of lattice and beam, such as equal beta functions and beta- 
tron tunes and beam sizes, in additional to maintaining the 
energy transparency condition: I+E + = I_E_ [Q]. The 
other one is to break some unnecessary symmetry, for ex- 
ample, betatron tunes. To make right choice, it is important 
to understand what are the consequences when symmetry 
is broken. For example, it is known that the violation of 
the energy transparency condition might cause a flip-flop 
of the colliding beams [ph. The main subject of this paper 
is to study the symmetry using the PEP-II as an example. 

First we simulate the beam-beam limit and the spectrum 
of coherent oscillation when the transparency conditions 
are preserved. Then we will study the spectrum and mo- 
tion of the coherent oscillation when the transparency con- 
ditions are largely violated. Finally we will make compari- 
son of the luminosity between the simulation and measure- 
ment. 

2 SYMMETRIC PARAMETERS 

Particle simulation is one of the important tools to study 
many aspects of the beam-beam interaction such as the 
beam-beam limit and the luminosity of colliders. Extend- 
ing the work [^, 0| of solving the Poisson equation, we 
reduce the region of mesh by assigning inhomogeneous 
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potential on the boundary [|lj]. The method allows us to 
choose much smaller region of mesh and therefore increase 
the resolution of the solver. The improved resolution makes 
more accurate the calculation of the dynamics in the core 
of the beams. 

In a typical simulation, we track 10240 macro particles 
inside an area of 8a x x 2Aa y with a rectangular mesh of 
256 x 256. For a beam aspect ratio of a x : a y — 32 : 1, 
we choose fifteen grids per a x and five grids per a y . This 
choice of simulation parameters makes about ten particles 
per cell on average within a region of 3a x x 3a y where the 
most of beam reside. It is adequate to compute the quanti- 
ties that are mostly determined by the core of the beam. 

The particles lost outside the meshed region are kept 
where they are lost and their contribution to the force is 
ignored afterward. The loss of the particles is closely mon- 
itored under different conditions since too much loss means 
that the simulated result is not reliable anymore. For a mesh 
size of 8173; x 24<jy, the loss of the particles is less than 1 % 
even at extremely high beam intensity. 

Due to the limitation of the computational speed on a 
computer workstation, only a single two-dimensional slice 
is used to represent a bunched beam. Therefore, all lon- 
gitudinal effects such as the hourglass effect and synch- 
betatron coupling are neglected in the simulations. 

At each beam intensity, we track the particles up to three 
damping time till the beams reach their equilibrium dis- 
tributions. Then the equilibrium distributions are used to 
compute the quantities like the beam-beam parameters. For 
extracting the power spectrum, we track additional 2048 
turns after the equilibrium and save the beam centroid ev- 
ery turn. 

2. 1 Parameters 

The PEP-II is an asymmetric e + e~ collider with two dif- 
ferent storage rings in a 2.2 kilometer tunnel at the Stanford 
Linear Accelerator Center. The positron beam is stored in 
the Low Energy Ring (LER); the electron beam in the High 
Energy Ring (HER). The two rings are vertically separated 
and brought into the collision at an interaction point (IP). 
In Tab. 1, we list a possible set of symmetric parameters 
for the PEP-II. 

It has been shown [[|] that even for an asymmetric col- 
lider, it still can be operated theoretically as if it is a sym- 
metrical collider if the transparency conditions: v+ = 

= Vy,P? = = %-,4 = = 

e~, r t + = r t , I+E+ = I-E-, are all satisfied as far as the 
beam-beam interaction is the concern. 

To preserve the energy transparency condition in the 
simulation, we vary the beam intensity with a step of 
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Table 1: Symmetric parameters for the PEP -11 
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10 10 and 5N~ = SN + E + /£__ starting from zero. 



2.2 Beam-Beam Parameters 

Given equilibrium distributions that are close enough to the 
Gaussian, we introduce the beam-beam parameters 



(1) 



where r e is the classical electron radius, 7 is the energy of 
the beam in unit of the rest energy, and N is total num- 
ber of the charge in the bunch. Here the superscript "+" 
denotes quantities corresponding to the positron beam and 
"— " quantities corresponding to the electron beam. 




Figure 1 : The beam-beam parameter as a function of single 
bunch current. The circles represent the positron beam and 
the cross represent the electron beam. The solid lines are 
fitted curves. 

At every the intensity of beams, we computed the beam 
sizes and the beam-beam parameters of the equilibrium 
beam distribution. The result is summarized in Fig. [I]. 
Clearly, the beam-beam parameter in the vertical plane is 
saturated at high beam intensity. That is consistent with 
many experimental observations [§]. Moreover, the depen- 
dence of the beam-beam parameter upon the single bunch 
current can be fitted rather well with two parameters, 



where £3,(00) is the beam-beam limit and a is the decay 
rate with respect to the beam current /. For this particular 
case, Cy(oo) = 0.0422. At I + = 1.26mA, which is the 
nominal value of the design single bunch current, £ y = 
0.025. That is about 15% less than the design value of 
the beam-beam parameter. At I + = 2.31mA, which is the 
single bunch current at the top of each filling last October, 
£ y = 0.033, which is less than £ a (oo). That indicates that 
there is still room to improve if the symmetric parameters 
could be implemented in the machines of the PEP-II. 

At a few other tunes that we have studied, for example 
v x = 0.2962 and v y = 0.2049, we found the similar phe- 
nomenon. The same is true in the horizontal plane. It is 
intriguing that such simple parameterization can be applied 
to the beam-beam parameter. 

2.3 Spectrum 

Using the fast Fourier transformation (FFT), we computed 
the power spectrum with the beam centroids which were 
recorded in 2048 consecutive turns after the equilibrium 
distributions were established. The several spectra with at 
different beam intensities are shown in Fig. |2|. There are 
two peaks clear seen in each spectrum. They are a and it 
modes of the coherent dipole oscillations. The tune shift of 
the 7r mode increases with respect to the beam intensity. 
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Figure 2: Power spectra of coherent dipole oscillation at 
different beam intensities. The dashed line represents the 
machine tune v x = 0.39. 

However, the tune shift of the a mode also increases 
though less than the n mode as the intensity increases. 
This phenomenon is beyond the capability of the linearized 
Vlasov theory since it predicts no tune shift for the a mode. 

2.4 The Yokoya Factor 

Studying the power spectrum of colliding beams is a pow- 
erful way to investigate and understand the beam-beam in- 
teraction. Historically, in symmetric colliders where two 
beams are identical, the tune shift of the coherent tt mode 
has provided many useful insights into the dynamics of the 
beam-beam interaction. It has been shown analytically that 



this tune shift is proportional to the beam-beam parameter 
£, namely 5u„ = A£ [g 0, [ll|, The coefficient A is 
between 1 and 2 depending on the beam distribution. For a 
self-consistent beam distribution JT^], 



A£ x ,5u yi7T = A(l - r)£ y 



(3) 



where A = 1.330 - 0.370r + 0.279r 2 , r = Oj,/^ + <? y ), 
and (jj; and a y are the horizontal and vertical beam size 
respectively. 

Experimentally, this relation has been observed in many 
different colliders Jol [l4| ]. The results of measurements 
are consistent with the calculation based on the Vlasov the- 
ory. In simulations fll5j, [T^ ] using the PIC this relation was 
also confirmed. Now, this well-established relation is often 
used to measure the beam-beam parameter or test a newly 
developed code. 

The shifts of the ir mode away from the machine tune 
are extracted from the spectra as shown Fig. ^ at different 
beam intensities. The shifts as a function of the beam-beam 
parameter, which is shown in Fig. [l], are summarized in 
Fig-! 




Figure 3: The tune shift of coherent n mode as a function 
of the beam-beam parameter. The left plot is for the hor- 
izontal plane and right plot is for the vertical plane. The 
circles represent the simulated tune shifts. The solid lines 
represent 8v v = A£. 

The predicted linear relation based on Eqn. is also plot- 
ted in the figure. One can see that the agreement between 
the theory and simulation is rather good even at very high 
beam-beam parameter. 

3 ASYMMETRIC PARAMETERS 

In October 2000, PEP-II has achieved its design luminosity 
of 3.0 x 10 33 cm~ 2 s -1 . The parameters of the lattices and 
beams when the design luminosity was reached are tabu- 
lated in Tab. 2. It is clear that many transparency con- 
ditions are violated. Among them, the betatron tunes are 
very different and well separated compared to the beam- 
beam parameter. 

During the last run, the ratio of the beam current I + : 2_ 
is about 2:1. As a result, the energy transparency condition 
I+E + = I-E- is also violated. With this set of parame- 
ters, the PEP-II has been operated at a region of asymmetry. 

To simulate the beam-beam effects under the nominal 
running condition of the PEP-II, we vary the beam inten- 
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Table 2: Operating parameters for PEP -11 

sity with a step of SN+ = 10 10 and SN~ = SN+ /2. Since 
the LER has longer damping time than the one in HER, 
we track the particles with three damping time of the LER 
to ensure that both beams reach their equilibrium distribu- 
tions. 

3. 1 Dipole Motion 

The Poincare maps of self-excited coherent dipole in 
the plane with four different beam intensities: N + = 
(1, 5, 9, 13) x 10 10 (from left to right) are shown in Fig. f| 



Figure 4: Poincare map of coherent dipole at different 
beam intensities in the horizontal plane. 

It is clear that there is maximum amplitude of the oscil- 
lation. The amplitudes are very small and within cr Xi j,/30 
when the beam intensities stay bellow the peak operating 
intensity N+ = 10.6 x 10 10 . 

In Fig. ^j, some structures of the resonance can be clearly 
seen. For instance, at N + = 1 x 10 10 , we see seven islands 
in the Poincare map for the electron beams. This seventh 
order resonance can be identified as lv~ = 172. For the 
positron beam, we see a triangle shape which is consistent 
with the third order resonance 3v x = 116 at N + = 5 x 
10 10 . 

It is worthwhile to note that the resonance structure dis- 
played in the figure is near a x /30, which is a factor of two 
smaller than the size of the grid (a x /l5). This does not nec- 
essarily mean that the resolution of resolving the dynamics 
of individual particle is less than the size of the grid. But 
for the collective motion, such as the oscillation of the co- 
herent dipole, the resolution can be smaller than the size 



of the grid. Because they are an average over the distribu- 
tion of the particles, the noises from the finite size of the 
mesh and the representation of beam distribution as a finite 
number of particles are much reduced. 

Below the peak operating intensity, the amplitude of the 
vertical oscillation is about a y / 30. But no higher order res- 
onances are identifiable. One of possible reasons is that the 
vertical grid size a y /5 is too large to resolve the structure 
of resonance within er y /30. 

Above the peak operating intensity, the amplitude of 
the vertical oscillation increases more than ten times and 
reaches half of the beam size for the electron beam. The 
oscillation acts coherently as a single particle. It is not clear 
why the dipole mode of the electron beam is excited to such 
large amplitude at the higher intensity. 

3.2 Spectrum 

In the vertical plane, the dipole mode at different beam in- 
tensities are plotted in Fig. |j. Unlike the symmetric col- 
liding beams, there are no visible it modes in the spectra. 
An asymmetric shape of the spectrum is clear visible es- 
pecially at the high intensity. Similar spectrum has been 
observed experimentally [17|. 
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Figure 5: The vertical power spectra at different beam in- 
tensities for the PEP-II. 

At N + = 13 X 10 10 , a single mode is coherently ex- 
cited in the electron beam. Correspondingly, the excited 
mode is shown as an ellipse in the Poincare map. And the 
positron beam is blowup vertically at the same time. As 
a result, the luminosity decreases. In order to check if this 
highly excited mode is the cause of the rapid increase in the 
vertical size of the positron beam, we eliminate the dipole 
oscillation every turn in simulation. But the peak luminos- 
ity remains the same. So we conclude, in this case, that 
the collective dipole motion is not the main reason for the 
beam-beam blowup. 

In the horizontal plane, the spectra are broader than the 
vertical ones largely because the resonances, which we 
have shown in the previous section. The spectra are not 
so asymmetric as the vertical ones. 



3.3 Tune Shift 

Similar to the symmetric case, the spectrum shifts as the in- 
tensity increases. The tune shifts as a function of the beam 
intensities is shown in Fig. |[ The center of the spectrum 
is the fitted result of the Lorentz spectrum. In general, the 
center does not coincide with the peak in the spectrum due 
to the asymmetric nature of the spectrum. Therefore, the 
tune shifts as plotted in the figure should be considered as 
the average values. The tune shifts saturate around 0.015 
in the both planes. In particular, the vertical shifts actually 
start to decrease near the peak operating intensity. Simi- 
lar behavior had been observed in the measurements of the 
power spectrum for the PEP-II. This behavior is certainly 
very different with the behavior of the beam-beam param- 
eters as simulated [jj]]. For the PEP-II, which is operated 
at very asymmetric parameters, we do not have a simple 
linear relation between the beam-beam parameter and the 
tune shift of the dipole spectrum. 




Figure 6: Tune shift as a function of bunch intensity. The 
circles represent the tune shift of the positron bunch. The 
crosses represent the tune shift of the electron bunch. The 
dashed lines represent the peak bunch intensity of the PEP- 
II operation. 



3.4 Luminosity 

To make a direct comparison between simulation and ex- 
perimental observation, we have recorded the luminosity 
during a period of four hours on October 1, 2000. The 
data are shown in Fig. ^[ Duration of each measurement 
was three minutes. The first and second plots in the fig- 
ure present the total decaying beam current of positron and 
electron beams respectively. The third plot shows the mea- 
sured and simulated luminosities at the same beam current 
displayed in the figure. The other parameters used in the 
simulation are the same as in Tab. 2. 

The agreement of the simulation and measurement was 
within 10%. Since the longitudinal effects of the beam- 
beam interaction are not yet included in the simulations, 
three-dimensional simulation could reduce the simulated 




Figure 7: Luminosity of a routine operation of PEP-II. The 
crosses represent measurement and the circles represent 
simulation. The number of bunches was 605. 

luminosity. For example, the hourglass effect should re- 
duce the simulated luminosity by 12% given a z = 1.3cm 
and {3* = 1.25cm. 

4 CONCLUSION 

When the transparency conditions are violated, especially 
the betatron tunes are well separated, the strong-strong sim- 
ulations show the spectra of coherent oscillation in the 
beam-beam interaction is very different from the spectra 
seen in the symmetric collider. In particular, there is no 
7r mode seen in the spectrum. The simple linear relation 
between the beam-beam parameter and the tune shift of 
the 7r mode is no longer existing. Given the operating pa- 
rameters of the PEP-II, we do not see any simple relation 
between the tune shift of the continuum spectrum and the 
beam-beam parameter. Therefore, the beam-beam parame- 
ter would not be estimated using the spectrum. 

The agreement of luminosity between the simulation and 
measurement is surprising and remarkable considering the 
simplicity of the two-dimensional model. In general, the 
three-dimensional effects such as the hourglass effects and 
the synch-betatron resonance could become very impor- 
tant. The code is being extended to include the bunch 
length and synchrotron oscillation. More simulation results 
will be directly compared to the controlled experiment. 
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